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' Abstract 

Reduced ordering based vector filters have proved successful in removing long-tailed noise from color images while 
preserving edges and fine image details. These filters commonly utilize variants of the Minkowski distance to order the color 
O |j vectors with the aim of distinguishing between noisy and noise-free vectors. In this paper, we review various alternative 
distance measures and evaluate their performance on a large and diverse set of images using several effectiveness and 
m : efficiency criteria. The results demonstrate that there are in fact strong alternatives to the popular Minkowski metrics. 

in ; 

i — 1 Introduction 
> 

Color images are often contaminated with noise, which is introduced during acquisition (sensor noise) and/or transmission 
^ (channel noise). Sensor noise is inherent to all electronic image sensors and is often modeled using Gaussian or Poisson 
(J distributions pQ. Channel noise is caused by various sources including man made phenomena (car ignition systems, industrial 
l—— 'machines in the vicinity of the receiver, switching transients in power lines, various unprotected switches, etc.) and natural 
I causes (lightning in the atmosphere, ice cracking in the Antarctic region, etc.) [2j- Channel noise is often characterized by 
impulsive sequences, which occur in the form of short-duration, high-energy spikes attaining large amplitudes with probability 
C — higher than that predicted by the Gaussian density model |3j . The introduction of such noise into an image not only lowers 
its perceptual quality, but also makes subsequent tasks such as edge detection and segmentation more difficult. Therefore, 
noise removal is often an essential preprocessing step in many color image processing applications. 
. ■ Numerous filters have been proposed for the removal of long-tailed, in particular impulsive, noise from color images 
[I] El El HI 16] . Among these, nonlinear vector filters have proved successful in removing the noise while preserving the edges 
and fine image details. The early approaches to nonlinear filtering of color images often involved the application of a scalar 
I filter to each color channel independently. However, since separate processing ignores the inherent correlation between the 
color channels, these methods often introduce color artifacts to which the human visual system is very sensitive. Therefore, 
. £^ vector filtering techniques that treat the color image as a vector field and process color pixels as vectors are more appropriate. 
An important class of nonlinear vector filters is the one based on robust order-statistics. These filters involve the ordering 
^ of a set of input vectors inside a predefined sliding window. The lowest ranked input vector in this window is often taken as 
■ - - the output vector. The ordering is commonly achieved using variants of the Minkowski distance. Recently, fuzzy similarity 
measures [71 [HI HI [TOJ [TTJ [T2J [T31 HH HU US] have been proposed as alternatives to the more traditional distance measures. 
Since it is often trivial to convert a similarity measure to a dissimilarity measure, throughout this paper we will use the term 
'distance measure' to refer to a function that quantifies the distance or similarity between two vectors. 

In this paper, we review various distance measures that can be used to order color vectors in reduced ordering based 
vector filters. We evaluate the performance of these measures on a large set of images that cover a variety of domains using 
several effectiveness and efficiency criteria. The rest of the paper is organized as follows. Section [2] presents an overview 
of order-statistics based vector filters and distance measures that can be utilized to order color vectors. Section [3] describes 
the noise model, filtering performance criteria, and the experimental setup. Finally, Section [4] gives the conclusions and the 
future work. 



2 Distance Measures for Reduced Ordering of Color Vectors 

Consider an MxN red-green-blue (RGB) input image X that represents a two-dimensional (2D) array of three-component vec- 
tors x(r, c) = [xi(r, c),X2(r, c),X3(r, c)] occupying the spatial location (r, c), with the row and column indices r = {I , . . . , M} 
and c = {1, . . . ,N}, respectively. In the pixel x(r, c), the Xk(r,c) values denote the red (k = 1), green (k = 2), and blue 
(k = 3) components. Most images are nonstationary in nature; therefore filters often operate on the assumption that the 
local image features can be extracted from a small image region called a sliding window. The size and shape of the window 
influence the properties and efficiency of the filtering operation and are therefore application dependent [T] . A square window 
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W(r,c) of size s/n x yfn pixels centered on x(r, c) is commonly used due to its versatility and good performance [I]. The 
window slides over the entire image X in a raster fashion and the procedure replaces the input vector x(r, c) with the output 
y(r, c) = F(W(r, c)) of a filter function F(-) that operates over the samples inside W(r, c). Repeating the procedure for each 
pair (r, c), with r = {1, . . . , M} and c — {1, . . . , N}, produces the output vectors y(r, c) of the M x N filtered image Y. For 
notational simplicity, the input vectors inside W(r, c) are re-indexed as a set, i.e. W(r, c) = {x 4 : i = 1, . . . , n} (sec Figure[T|), 
as commonly seen in the related literature [TJ [5J [31 [U 03 [B] . In this notation, the center pixel in W is given by X( n+1 )/ 2 an d 
in the vector Xj = [xn, Xn, 2^3] with components x^, the i 6 {1, . . . , n} and k £ {1, 2, 3} indices denote the position of the 
vector inside the window and the color channel, respectively. 

Xj x 2 x 3 

X 4 X 5 X 6 
_ X 7 X 8 X 9_ 

Figure 1: Indexing convention inside a 3 x 3 window 

Order-statistics based vector filters operate by ordering (ranking) the vectors inside the sliding window. The purpose 
of this ordering is to distinguish between the noisy and noise-free vectors. However, in contrast to scalar data, there is no 
universally accepted method for ordering multivariate data. Widely known multivariate ordering methods include |17j : 

• Marginal ordering (M-ordering): The vectors are ordered in each component independently This scheme produces a 
set of ordered output vectors that is usually not the same as the set of input vectors and consequently results in color 
artifacts when applied to multichannel image data. 

• Conditional ordering (C-ordering): The vectors are ordered based on the marginal ordering of one of the components. 
This scheme disregards the vectorial nature of the multichannel image data. 

• Partial ordering (P-ordering): The vectors are partitioned into smaller groups that are then ordered. Despite its 
theoretical appeal, this scheme is computationally demanding |18j . 

• Reduced (aggregate) ordering (R-ordering): The vectors are first reduced to scalar representatives using a suitable 
distance measure. The ordering of these scalars is then taken as the ordering of the corresponding vectors. This is the 
most common ordering scheme in the literature [TJ [21 El 01 [6] . 

To order the vectors inside the sliding window W(r, c) — {x^ : i = 1, . . . , n}, R-ordering based vector filters first calculate 
the aggregate (cumulative) distance Di for each vector x^: 

n 

A = 5>(x,:,x,) (1) 

3=1 

where d(xj,Xj) denotes a suitable distance measure that quantifies the dissimilarity between vectors Xj and Xj. The ordered 
sequence of scalars < D( 2 ) < . . . < -D(„) implies the same ordering in the corresponding vectors x^) < X( 2 ) < ■ ■ ■ < X( n ). 
Many filters define the lowest ranked input vector (lowest order-statistic), i.e. x^), as the output vector. This is because 
vectors with lower ranks are typical (representative) vectors in their neighborhoods, whereas those with higher ranks are 
often outliers (noisy vectors) [19] . More explicitly, the output vector at location (r, c) is given by: 

y(r, c) = argmin Di (2) 

XigW(r,c) 

Note that if a similarity measure s(xj,Xj) is used instead of d(xj,Xj) in (p}, the argmin operator should be replaced with 
argmax. 

The vector median filter (VMF) [50], which is the most well-known R-ordering based vector filter, utilizes the Minkowski 
/ 3 \1/p 

distance dp(xi,Xj) = [J2k=i \ Xik ~ x jk\ P ) as its ordering criterion. Three special cases of d p are of particular interest: 
d\ (city-block distance), d 2 (Euclidean distance), and doo (chessboard distance). Given the general Minkowski form, d\ and 
di can be defined in a straightforward fashion, while d^ is defined as doo(xi, Xj) = maxi<fe<3 \xik — Xj%\. Among these three 
forms, di seems to be preferred in the color image filtering literature due to its isotropic (rotation invariant) nature and good 
performance [21]. The squared form of di, i.e. d|, can also be used as a computationally efficient alternative [22] . 

Table Q] lists some of the traditional distance measures that can be used for the reduced ordering of color vectors. Note 
that the Bray-Curtis, Canberra, Soergel, and Ware-Hedges measures are normalized forms of the city-block distance, whereas 
the Divergence Coefficient measure is a normalized form of the Euclidean distance. Some of the measures given in Table [TJ 
have already been used in the design of various R-ordering based vector filters. For example, the Cosine distance was used 
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Table 1: Distance measures for R-ordcring of color vectors (||x||: Euclidean norm of x) 



Bray-Curtis Distance [26] 



Canberra Distance [27] 



Chord Distance [28] 



Cosine Distance l23l 



Divergence Coefficient [29] 



Goude Distance 



Soergel Distance 



Ware-Hedges Distance [32] 



E \x ik -x jk \ 

<^bray (Xj , Xj ) 



E Xik+Xjk 
fc=l 



^Canberra (Xi , X^ ) — a; ifc +a; 
fc=l 



2* ik j k | 

3k 



^chord (Xi , Xj ) — | ("\/ ^z/c y/^jk 

\k=l 



E x ik X jk 

^cosine ( x i , Xj ) = 9 = ai'CCOS j -jj^j y || 



rfdivcrgcncc(Xi, Xj) = ( S x ik +x jk 



|x,|| 2 + ||x 3 || 2 -2 J2 x ik x jh 

^goudc (Xj , Xj ) 



2 + ||x J || 2 +2 E XifcXjfc 
fc=l 



E l^ifc-^jt 

^soergel (Xj , Xj ) 



E max(i it ,i jt ) 
fc=l 



dv/are (Xj , Xj ) — 



fc=l 



\xjk—Xjk\ 

max(x ik ,x jk ) 



in the design of the basic vector directional filter (BVDF) [53]. The Goude distance was used in the content-based rank 
filter [21]. The rationale behind this measure is that the similarity between two vectors x» and Xj can be expressed as the 
ratio of some function of what they share (commonality) to what they comprise (totality). In particular, the numerator and 
denominator correspond to the vector difference and the vector sum, respectively. 
In [9], the authors introduced a fuzzy magnitude similarity measure given by: 



s fms ( x i ! x j ) — J 



fc=i 



mm(xjk,Xj k ) + K 
max(xik,Xjk) + K 



(3) 



where K is a parameter of the measure. 

The analogous fuzzy directional similarity measure between vectors Xj and Xj can be obtained by first normalizing the 
vectors to unit length, i.e. x, = Xj/||xj|| and xj = Xj/||xj||, and then calculating the fuzzy magnitude similarity between 
them, i.e. sjy s (xi,Xj) = s|^ s (x,,Xj) [11]. Furthermore, since both fuzzy magnitude similarity and fuzzy directional similarity 
values fall into the (0, 1] interval, they can easily be combined to obtain a fuzzy magnitude-directional similarity measure as 
follows: sg£'(xi,Xj) = sj^fa,Xj) • s$£(xi,Xj) [TQ. 

In |10j . the authors introduced a combined fuzzy similarity measure given by: 



S cfs l Xi ' X J/ 



c 



C + d 2 (xi, Xj) t + max (\n — Tj\, \a - Cj 



(4) 



where (r,,Ci) and {rj,Cj) are the spatial coordinates of the vectors Xj and Xj, respectively, and C and t are the parameters 



of the measure. It can be seen that s cf '* combines color similarity (first term) and spatial proximity (second term). The 
maximization of the aggregate combined fuzzy similarity in a particular neighborhood is equivalent to finding the vector that 
is centrally located and is a good representative of the neighborhood in terms of chromatic content. 

The ordering criterion used in the directional-distance filter (DDF) [23] combines the Minkowski and Cosine distances. In 
particular, the aggregate distance associated with pixel x, is given by: 



D{ I ^ ^ dp (x^ , Xj ) I I ^ d cos [ nc (x^ , Xj ) 



(5) 



Note that the idea of combining directional and magnitude processing also underlies the design of Sf m( js- However, while Sf m d s 
is a proper similarity measure, i.e. it quantifies the similarity of two vectors, the ordering criterion given in ([5]) is only defined 
in a particular neighborhood. 

A function / : S x 5* — > R is called a metric on a set S if it satisfies the following conditions for all x,y,z £ S [55] : 
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Table 2: Operation counts for the distance measures 



Measure 


ABS 


COMP 


ADD/SUB 


MULT 


DIV 


SQRT 


ARCCOS 




3 




5 
















5 


3 




1 




a 2 






5 


3 








rf 


3 


2 


3 










"bray 


3 




10 




1 






^Canberra 


3 




s 

o 




3 






'-'CIS 


2 


x 


9 


4 


2 


1 




"'chord 






5 


3 




7 




rl ■ 

"'cosine 






6 


10 


1 


1 


1 


"-divergence 






s 

O 


3 


3 


I 




s fds 




3 


1 


« 

o 


9 


2 








3 


6 


2 


3 






3fmds 




6 


16 


11 


12 


2 




^goudc 






9 


10 


1 


1 




^soergel 


3 


3 


7 




1 






^ware 


3 


3 


5 




3 







1. f(x,y) > (non-negativity) 

2. /(a;, y) = if and only if x = y (self-identity) 
3- f(x,y) = f{y,x) (symmetry) 

4. f{x, y) < f(x, z) + f(y, z) (triangle inequality) 

Among the distance measures presented in this section, the Squared Euclidean, Bray-Curtis, and Goude measures are not 
metrics on K since they violate the triangle inequality. The Cosine distance is a semi (pseudo) metric on K because it satisfies 
all of the conditions except for the self-identity condition. On the other hand, it can be shown that [9j [TUJ [TT| the fuzzy 
similarity measures Sf ms , Sfd s , Sf m( js! an d s c f s are stationary F-bounded fuzzy metrics in George and Veeramani's sense |34j . 

The operation counts for each distance measure are given in Tabled (ABS: absolute value, COMP: comparison, ADD: 
addition, SUB: subtraction, MULT: multiplication, DIV: division, SQRT: square root, ARCCOS: inverse cosine). Note that 
these operations are listed in ascending order according to their computational costs on a typical processor architecture, 
i.e. cost(ABS) < cost(COMP) < ...< cost(SQRT) < cost(ARCCOS). This may not hold on certain architectures, but it 
nevertheless permits cost comparisons between the measures to a certain extent. The actual cost of each operation is highly 
dependent on the implementation platform (the processor and the compiler). 

It should be noted that, even if it is mathematically simple, the use of a particular distance measure in an R-ordcring 
based vector filter can be computationally demanding. This is because during the filtering procedure the aggregate distance 
value |T]) is calculated for every vector Xj in each neighborhood W(r,c). This means that the determination of the output 
vector y(r, c) requires n(n — l)/2 distance evaluations. For anMxJY image this amounts to MNn(n— l)/2 evaluations. 
As an example, filtering a 512 x 512 image using a 3 x 3 window requires over 9.4 million evaluations. Fortunately, the 
computation of certain distance measures can be accelerated using table lookups (see SJ3]). 



3 Experimental Results and Discussion 



In order to evaluate the performance of the presented distance measures, a set of 100 high quality RGB images was collected 
from the Internet. The set included images of people, animals, plants, buildings, aerial maps, man-made objects, natural 
scenery, paintings, sketches, as well as scientific, biomedical, synthetic, and test images commonly used in the color image 
processing literature. The corruption in the images was simulated by the widely used correlated impulsive noise model 35 : 



x = < 



o with probability 1 — , 

{r\, 02, 03} with probability ipi ■ ip, 

{01, r%, 03} with probability tp% ■ tp, 

{01, 02, r^} with probability (^3 • ip, 

{ri,r 2 ,r 3 } with probability (1 - (tpi + ip 2 + ^3)) • <p 



(6) 



where o = {01,02,03} and x = {x\, x 2 , £3} represent the original and noisy color vectors, respectively, r = {ri,r2,r3} is 
a random vector that represents the impulsive noise, ip is the sample corruption probability, and ipi, ip 2 , and ip 3 are the 
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corruption probabilities for the red, green, and blue channels, respectively. In the experiments, the channel corruption 
probabilities were set to 0.25. 

Filtering performance was evaluated by the following effectiveness criteria [BJ: 

M N 

1. Mean Absolute Error: MAE (X, Y) = 3^ E E d i ( x <7, c), y(r, c)) 

r— 1 c— 1 

where X and Y denote respectively the M x N original and filtered images in the RGB color space. MAE measures 
the detail preservation capability of a filter. 

M N 

2. Mean Squared Error: MSE (X, Y) = 3^ E E <% (x(r, c), y(r, c)) 

r— 1 c— 1 

MSE measures the noise suppression capability of a filter. 

M N 

£ £ d2( X Lab (r,c),y Lsb (r,c)) 

3. Normalized Color Difference: NCD(X,Y) = r=lc=1 M N 

£ £ 11X^(7^)11 

r=l c=l 

where x Lab (r, c) and y Lab (r, c) denote the CIEL*a*b* coordinates of the pixel (r, c) in the original and filtered images, 
respectively. NCD measures the color preservation capability of a filter. 

Computational efficiency was measured by CPU time in milliseconds. All of the programs were implemented in the C 
language, compiled with the gec 4.2.4 compiler, and executed on an Intcl®Corc™2 Quad Q6700 2.66 GHz machine. The 
time measurements were averaged over 10 identical runs. Although relative efficiency judgments can be made using the 
information given in Table [5J the actual execution times were measured for two main reasons. First, the cost of the basic 
operations is not uniform (for example, ARCCOS is computationally very expensive when compared to the others). Second, 
the use of lookup tables (LUTs) can significantly accelerate the computation of certain distance measures. 

Eighteen distance measures were evaluated in the experiments. These included the four variants of the Minkowski distance 
(di, g?2, doc, and c?|), the eight measures listed in Table [1] and the four fuzzy similarity measures discussed in Sj2] Each of 
these measures was plugged into (TTJ) and the performance of the resulting filter was evaluated using the abovementioned 
criteria. A square window of size 3 x 3 (n = 9) is used in all of the filters. Although the ordering criterion used in DDF ([5]) is 
not a distance measure, because of its widespread use in the literature, it was included in the experiments. Furthermore, the 
identity filter (indicated by the label 'none'), which performs no filtering was included as a baseline. The K parameters of 
Sf ms and Sfds were set to 1024 and 4 [11], respectively. The C and t parameters of s c fs were set to 150 and 4 [10], respectively. 

Implementation details for the distance measures are given below (0 = [0,0,0]): 

> d\: Instead of calculating n (n — l)/2 distance values in each neighborhood W(r, c), first the mean vector x = l/<T-E™=i Xi 
is calculated and then the output vector is determined to be the input vector that has the minimum d\ distance to the 
mean vector, i.e. y(r, c) = argmin d|( x i! x ) |17| . 

Xi£W(r,c) 

> dbrayj ^goudc, d socrgc i: If Xj = Xj = 0, the distance between the two vectors is taken as 0. 

> d can bcrra, ^divergence, dware ■ If ^ifc = Xjk = for any component k € {1, 2, 3}, the contribution of that component to the 
distance between Xj and Xj is ignored. The value of the term in the summation is precomputed and stored in a 2D 
LUT of size 256 x 256. 

> s c f s : The spatial proximity term is precomputed and stored in a 2D LUT of size n x n. 

> dchord: The value of each term in the summation is precomputed and stored in a ID LUT of size 256. 

> dcosine: If either x, = or Xj = 0, the distance between the two vectors is taken as tt (maximum possible angular 
distance). If Xj = Xj = 0, the distance between them is taken as (minimum possible angular distance). 

> Sfd s : If either x, = or Xj = 0, the similarity between the two vectors is taken as (minimum possible similarity). If 
x, = Xj = 0, the similarity between them is taken as 1 (maximum possible similarity). 

> Sf ms : The value of the term in the product is precomputed and stored in a 2D LUT of size 256 x 256. 

Table [3] shows the average effectiveness ranking of the distance measures over the entire image set. Note that the ranks 
start from and the smaller the rank value, the better the corresponding distance measure. For example, with respect to 
the MAE criterion, the Sf m( js measure has an average rank of 2.00 at 20% noise level. In other words, with respect to detail 
preservation, the R-ordering based vector filter that utilizes Sf m( js ranks, on the average, in top 3 among 18 filters. The last 
column shows the overall mean effectiveness rank for each distance measure. It can be seen that, on the average, s c f s , Sf m d s , 
and Sf ms perform the best. Interestingly, these are all fuzzy metrics, which indicates that fuzzy logic might be better suited 
for outlier detection when compared to traditional distance measures. Among these metrics, s c fs outperforms the other two 
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by a large margin as evidenced by its impressive 0.33 average rank. The success of s c f s is most likely due to its consideration 
of spatial proximity, an idea that seems to be overlooked in the literature. It should also be noted that the top two filters 
are both hybrid in nature: s c fs combines magnitude similarity with spatial proximity, whereas Sf m ds combines magnitude 
similarity with directional similarity. In addition to the fuzzy metrics, d\ and divergence also perform better than c^, which 
is the most commonly used distance measure in the literature [1] [U |3l |U [5j |6] . Despite the fact that it differs from c?2 only 
by its lack of a square root operation, d\ performs among the worst. This is not surprising given the fact that the R-ordering 
based vector filter that utilizes d\ essentially approximates the average filter [55], which is known to blur the fine details, 
^cosine and Sfds also perform poorly, which shows that directional information alone is not sufficient to determine the most 
representative vector in a neighborhood. Interestingly, the identity filter ('none'), which leaves the noisy image unchanged 
preserves the fine details better than Sfd s , d%, and d C osine at 10% noise level. 

It should be noted that the performance of the fuzzy metrics can be enhanced by tuning their parameters adaptively 
according to the image characteristics and the level of noise, see, for example, |10j . 

Figure [2J shows sample filtering results for a closc-up part of the Peppers image contaminated with 10% noise. It can be 
seen that d\ blurs the fine details severely, (^divergence, £^2, d\, and Sf ms are fairly close in terms of detail preservation and noise 
removal performance. On the other hand, s c f s performs significantly better than the others, producing an output image that 
closely resembles the original image. 

Figure |3] shows sample filtering results for a closc-up part of the Lenna image contaminated with 20% noise. It can be 
seen that d C osine not only fails to remove a significant amount of noise, but also smears the fine details. Even though d gc ,ude, 
d2, Sfmds, and d\ suppress the noise well, this comes at the expense of the blurring of the fine details, e.g. the boa fur and 
the eye lashes. In contrast, s c f s achieves a remarkable balance between detail preservation and noise removal. 

Figure [4] shows sample filtering results for a close-up part of the Baboon image contaminated with 30% noise. This image 
presents a challenging case not only because of its high level of noise, but also due to its complex structural content, e.g. 
the whiskers. It can be seen that d c hord and c?2 fail to preserve the fine details, producing output images with blurry and/or 
broken whiskers. d W arc, dsoergel, d±, and Sf ms preserve the details better, while leaving many noisy pixels intact. As before, 
s c f s outperforms the others in terms of both detail preservation and noise removal. 

Table [4] shows the average efficiency ranking of the distance measures. Obviously, the identity filter consistently ranks 
first since it involves no actual filtering. As expected, d\ is the most efficient distance measure since, in contrast to the other 
measures, it does not involve n(n — l)/2 distance evaluations in each neighborhood. d\ is the second most efficient measure 
as it involves only the simplest mathematical operations, i.e. absolute value and addition/subtraction. The computation of 
d\ can be further accelerated using the method described in [3B|. d waro , Sf ms , and d C anberra are also among the most efficient 
since they can be calculated rapidly using LUTs. Except for Sf ms , the fuzzy metrics are not particularly efficient. Among 
these, s c f s and Sf m d s partially benefit from the use of LUTs. However, the calculation of the color similarity term in s c f s and 
the directional similarity term in Sf m d s increases the computational requirements for these filters. Finally, d C osine and ddf are 
the least efficient of all, which is due to their use of the computationally expensive ARCCOS operation. 

Table [5] shows the average computational time requirements of the distance measures on the 512 x 512 Peppers image 
over 100 identical runs. The last two columns give the actual (measured in seconds) and relative (indicated in units of t) 
computational times for each distance measure. It can be seen that the actual time differences among many of the measures 
are rather negligible. However, d cos ine and ddf arc still relatively time consuming when compared to the others. 

4 Conclusions and Future Work 

In this paper, wc reviewed various distance measures that can be used to order color vectors in R-ordcring based vector 
filters. We evaluated the performance of these measures on a large and diverse set of images using several effectiveness and 
efficiency criteria. The results demonstrated that the fuzzy similarity metrics s c fs, Sfmds, and Sf ms significantly outperform 
the commonly used d\ and di metrics in terms of detail preservation, noise removal, and color preservation. Furthermore, 
these fuzzy metrics are parametrized, which permits the tuning of their parameters according to the image characteristics 
and the level of noise. 

The use of the presented distance measures is not limited to the basic R-ordering based vector filters such as VMF, 
BVDF, and DDF. Each of these basic filters can in turn be used as a back-end in the implementation of more advanced 
switching filters, i.e. those filters that selectively remove the noise by utilizing impulse detectors [TJISJIS]. However, a distance 
measure that performs well in a non-switching noise removal filter may not necessarily perform as well when combined with 
an (imperfect) impulse detector. In other words, the complex interaction between the front-end (impulse detector) and back- 
end (noise removal filter) in a switching filter will likely to have a non-negligible influence on the ultimate performance of a 
distance measure. This issue will be investigated in a future study. 
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Table 3: Average effectiveness ranking of the distance measures 



Measure 


10% 


MAE 
20% 


30% 


10% 


MSE 
20% 


30% 


10% 


NCD 
20% 


30% 


Mean 


S c fs 


0.03 


0.00 


0.00 


0.00 


0.25 


2.67 


0.00 


0.00 


0.04 


0.33 


^fmds 


2.44 


2.00 


2.23 


3.42 


2.42 


2.01 


1.47 


1.21 


1.30 


2.06 


^fms 


2.20 


1.65 


1.30 


2.07 


1.79 


1.62 


6.49 


3.21 


2.53 


2.54 


di 


3.45 


2.86 


2.63 


1.99 


2.42 


2.61 


8.13 


4.90 


4.06 


3.67 


^divergence 


6.23 


6.08 


4.85 


4.67 


3.74 


2.94 


3.38 


3.25 


3.21 


4.26 


d 2 


7.61 


7.98 


7.24 


5.97 


5.01 


4.18 


4.12 


5.46 


5.34 


5.88 


^chord 


9.78 


9.69 


8.78 


7.86 


6.58 


5.53 


6.92 


7.66 


6.76 


7.73 


^Canberra 


6.10 


5.24 


6.03 


6.48 


8.35 


10.43 


10.32 


8.97 


10.46 


8.04 


ddf 


9.53 


9.99 


9.13 


11.00 


9.05 


7.28 


4.95 


5.85 


5.84 


8.07 


^bray 


6.03 


6.40 


8.26 


6.32 


9.19 


11.24 


10.16 


10.05 


11.14 


8.75 


^goude 


9.30 


11.08 


11.61 


8.66 


9.43 


8.94 


5.84 


8.85 


8.99 


9.19 


^soergel 


9.02 


7.64 


7.21 


10.25 


9.82 


10.46 


10.90 


8.89 


8.80 


9.22 


doc 


12.19 


12.59 


12.53 


11.11 


11.16 


9.25 


6.15 


10.81 


10.01 


10.64 


^warc 


10.49 


8.35 


9.39 


11.67 


12.36 


13.51 


12.95 


12.24 


13.12 


11.56 


Sfds 


14.75 


14.12 


14.09 


14.90 


14.73 


14.72 


13.93 


13.81 


13.65 


14.30 


d\ 


15.90 


15.73 


15.64 


14.11 


14.03 


12.76 


15.83 


15.88 


15.80 


15.08 


^cosine 


15.45 


15.12 


15.16 


15.52 


15.67 


15.85 


14.73 


14.96 


14.95 


15.27 


none 


12.50 


16.48 


16.92 


17.00 


17.00 


17.00 


16.73 


17.00 


17.00 


16.40 



Table 4: Average efficiency ranking of the distance measures 



Measure 


Time 
10% 20% 30% 


Mean 


none 


0.00 0.00 0.00 


0.00 


di 


1.00 1.00 1.00 


1.00 


di 


2.00 2.00 2.00 


2.00 


^ware 


3.48 3.42 3.47 


3.46 


£fms 


3.52 3.58 3.54 


3.55 


^Canberra 


5.00 5.00 4.99 


5.00 


doc 


6.00 6.01 6.00 


6.00 


^bray 


7.00 6.99 7.00 


7.00 


cfchord 


8.61 8.30 8.09 


8.33 


^soergel 


8.41 8.71 8.92 


8.68 


d 2 


10.02 10.01 10.08 


10.04 


^divergence 


10.98 11.00 10.93 


10.97 


°cfs 


12.06 12.05 12.06 


12.06 


^goudc 


12.96 12.97 12.96 


12.96 


Sfds 


13.97 13.97 13.97 


13.97 


Sfmds 


14.99 14.99 14.99 


14.99 


^cosine 


16.00 16.00 16.00 


16.00 


ddf 


17.00 17.00 17.00 


17.00 
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Table 5: Time requirements of the distance measures on the Peppers image 



Measure 


Actual Time 


Relative Time 




0.000 


0t 


d{ 


0.063 


It 




0.153 


2t 


3fms 


181 


3t 


d 

"ware 


0.182 


3t 


d u 

"Canberra 


183 

\J . -LOU 


3t 


"oo 


0.240 


4t 


"bray 


0.256 


4t 


d i 

"socrgcl 


0.346 


6t 


"chord 


356 


fiT 


"2 


366 


fiT 
U 1 


"'divergence 


0.396 


6t 


Sets 


508 


8t 


^goudc 


0.533 


8t 


Sfds 


0.623 


10t 


^fmds 


0.694 


llT 


^cosine 


1.424 


23t 


ddf 


1.684 


27t 



and presentation of this paper. 
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Figure 4: Filtering results for the Baboon image corrupted with 30% noise 
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